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Abstract 

The coefficient Cgw appearing in the Sheikholeslami-Wohlert improved action is 
computed to one loop perturbation theory for improved gluon actions including 
six-link loops. The 0(a) improvement coefficients for the dimension three isovector 
composite operators bilinear in the quark fields are also computed to one loop order 
of perturbation theory with degenerate non-vanishing quark masses. 
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1 Introduction 


Despite the fact that the quenched approximation is a drastic modihcation of the theory, 
it turns out that the low lying hadronic spectrum obtained in this approximation is 
quantitatively very close to the experimentally observed one. Only recently, with the 
progress of lattice QCD simulations, have the precise results from CP-PACS |I|, || shown 
a systematic deviation not only in the meson sector but also in the baryon masses. It is 
thus now clear that one has to incorporate the effect of dynamical quarks to have a chance 
to obtain the correct hadron spectra in the continuum limit. On the other hand, full QCD 
simulations are extremely computer time consuming compared to those of quenched QCD. 
Although the fastest computers now available have peak speeds of a few hundred Chops, 
it is still difficult to reduce the lattice spacing less than 0.1 fm while keeping the physical 
volume larger than 2 fm with reasonable statistics. Due to the large scaling violation 
of the combination of the plaquette gluon action and the Wilson quark action there is 
no hope to obtain accurate and reliable hadron spectra in the continuum limit by the 
extrapolation from such lattice spacings. 

One of the possible ways to overcome this difficulty is, of course, the use of improved 
actions to reduce the scaling violation. Among the several proposals, the least ambitious 
but most practical one known at present is the on-shell 0(a) improvement program, where 
only one new term (clover term) has to be added to the quark action (clover action) |^. 
In this program the scaling violation is reduced to O(a^) with the appropriate value of 
the coefficient of the clover term (clover coefficient). Although no modihcation for the 
plaquette gluon action is necessary in the 0(a) improvement program, one can still add 
one or more terms to the plaquette gluon action without signihcantly increasing CPU 
time, in order to reduce non-negligible O(a^) errors present at lattice spacings equal to or 
larger than 0.1 fm. Recently such a combination (the so-called renormalization group (RG) 
improved gluon action [Q and the clover action) has been compared to other combinations 
at a lattice spacing of ~ 0.2 fm and it has been shown that this combination has some 
better properties than others]^]: the rotational symmetry of the static quark potential is 
well restored and the scaling violation in hadron spectra is reasonably small. 

Because of this result, the CP-PACS collaboration decided to use the RG improved 
gluon action and the clover action for their production run of full QCD simulations 
whose lattice spacings range from 0.2 fm to 0.1 fm. Unfortunately at the time they 
started the production runs the value of clover coefficient Csw for the RG improved gluon 
action was known neither non-perturbatively nor at 1-loop order of perturbation theory. 
Instead they used a “perturbative mean held” value P Csw = (1 — 0.8412/3“^)“^^^, where 
1 — 0.8412/3“^ is the one-loop value for the plaquette. Although this value might be a 
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dominant contribution at 1-loop order, it is of course necessary to know Csw fully at 1-loop 
order to be able to estimate how large the errors in this approximation are. To obtain 0(a) 
improvement for correlation functions of composite operators one also has to construct 
0(a) improved operators [0. For the case of quark bilinear operators of dimension 3 (see 
subsection 2.4) this involves computing mixing coefficients (cx, X = A, V, T) as well as 
mass dependent correction factors (1 -|- amby, Y = S, P, A, V, T, m). It would of course 
also be desirable to push 0(a) improvement to the non-perturbative level as it has been 
done by the Alpha Collaboration for the plaquette gluon action. 

In this paper we have calculated Csw, cx, and by at 1-loop order of perturbation theory 
for gluon actions including six-link loops and the clover quark action. Our computational 
method uses the Schrodinger functional 0. Our results for the coefficients cx, and by 
are in complete agreement with the results of the recent computation of Taniguchi and 
Ukawa who used a completely different method. 


2 Definitions 


As mentioned in the introduction we are here concerned with constructing 0(a) on-shell 
improved actions for QCD which are to be used in large scale numerical simulations. 
Starting with the original Wilson action S'w for the quarks, Sheikholeslami and Wohlert 
0 have shown that 0(a) on-shell improvement can be achieved by adding just one extra 
term oc Fsw to the action. Thus we will be considering total actions of the form 

S[U,^,i;]=S[U] + Sy,[U,iP,'iP] + cMSsw[U,^,i^], ( 2 . 1 ) 


with Nf >2 flavors. The part of the action involving the quarks will be specihed in more 
detail below. At this point we only stress that the factor Csw multiplying the SW term 
is a function of the bare coupling which depends on the particular form of the gauge 
action 5[f/] chosen. The tree-level value is the same for all gauge actions Cs„(0) = 1 0 [] 
and so Cs^ has a perturbative expansion of the form 


Csw(fi'o) — 1 + C^wS'o + • • • 


( 2 . 2 ) 


For the Wilson plaquette action the 1-loop coefficient was computed a long time ago by 
Wohlert (and checked in ref. [|T^]), and recently non-perturbatively for the quenched 
theory |]^| and full QCD with Af = 2 flavors [|I4| . 


^Here we take Wilson’s r-coefHcient fixed to be 1, otherwise Csw(O) = r 
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In this paper we will compute for some gauge actions belonging to a general class 
containing loops up to length 6 

■Sl'^l = 4i;c<(So) E-C(C) (2.3) 

i=o c&Si 

where the Si denote sets of elementary loops C on the lattice as given in £g. |l|. Loops C 




( 0 ) 


( 1 ) 




( 2 ) 


(3) 


Figure 1: Elementary loops in Si, i = 0,1, 2, 3 

that differ by orientation only are considered equal. Furthermore 

C{C) = ReTr[/- t/(C)], 

^ Our notation differs from refs. by an interchange of C 2 with C 3 


(2.4) 
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U{C) being the ordered product of the link variables U^{x) along C. The coefficients Ci{gQ) 
are regular at = 0 and are normalized such that 

co(fi'o) + 8 ci((7o) + 16c2((7o) + 8 c 3 (fi'o) = 1- (2-5) 


This class of gauge actions was hrst considered by Wilson in the framework of 
block spin renormalization and includes in particular the action proposed by Iwasaki 0. 
Among this class are also the on-shell O(a^) improved pure gauge actions |^. However 
we stress that here we do not have the ambition to achieve on-shell O(a^) improvement 
for full QCD. The latter ambitious program would involve the addition of various other 
terms to the action including in particular 4-fermion operators (of dimension 6 ) which 
would be very awkward for numerical simulations. 

In our explicit computations we have only considered actions with C 2 = 0. Apart from 
the Wilson action (ci = C 3 = 0), we have studied 4 actions 


LW : ci(0) = 

- 1 / 12 , 

C3(0) = 0 

|161. 

(2.6) 

RGl : Cl = 

-0.331, 

C3 = 0 i 


(2.7) 

RG2 : Cl = 

-0.27, 

C3 = —0.04 

i, 

(2.8) 

RG3 : Cl = 

-0.252, 

C 3 = -0.17 

0. 

(2.9) 


Note for the LW action the Cj are functions of the bare coupling, whereas for the RG 
actions the q are constants. 

In ref. [^] it was shown that demanding positivity of the action in the limit 0 

restricts the tree level coefficients to some convex domain P. A sufficient condition for 
the coefficients to lie in P is 


where 

Ci = |(ci(0) - 

Note that all the actions listed above satisfy this condition except RG3. 


80. 

—C 3 > 0 

3 

(2.10) 

(0)1). 

(2.11) 


2.1 The Schrodinger Functional 

The work of the Alpha Gollaboration has demonstrated that a convenient framework 
in order to compute improvement coefficients and renormalization constants (and of 
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course more importantly many quantities of direct physical relevance) is provided by 
the Schrodinger functional (SF) ^ Here one considers the theory dehned on hy- 

percubic lattices of volume xT with cylindrical geometry, i.e. periodic-type boundary 
conditions in the spatial directions and Dirichlet boundary conditions in the “time” di¬ 
rection, for both gluon and quark helds From now on we set the lattice spacing 

a = 1 in most formulae. 

H, 0. 


For the case of the Wilson action the SF was described in detail in refs. 


The extension of the SF considerations to gauge actions belonging to the class above 
was hrst considered by Klassen m. We found it conceptually more convenient to set 
things up slightly differently. In particular in our construction the dynamical variables 
to be integrated over are chosen independently of the form of the action: they are the 
spatial link variables and quark helds Uk{x),'ip{x) with times xq = 1,.., T — 1 and the link 
variables Uq{x) with times Xq = 0, ..,T — 1 (i.e. inside the cylinder). Dirichlet boundary 
conditions are imposed on the helds Uk{x),'ip{x) at the boundaries Xq = 0 and Xq = T. 

What remains to be done is to specify the action in more detail. 


2.2 The pure gauge action S[U] 

Taking account of the fact that we are working with the Schrodinger functional geometry, 
we modify the gauge action above to read 

= ( 2 . 12 ) 

9o i=oc&Si 

with weights Wi{C,gQ) which may diher from Ci{gQ) for loops C near the boundaries. 
First we must specify the elements of the classes Si more precisely. These consist of all 
loops of the given shape which can be drawn on the cylindrical lattice that involve only 
the “dynamical links” in the sense specihed above, and spatial links on the boundaries at 
xo = 0 and xq = T. In particular rectangles protruding out of the cylinder are not included 
and hence we do not have to specify further boundary conditions for link variables outside 
the cylinder. 

The most important point now is that to achieve 0(a) improvement for all (on-shell) 
Green functions, a careful choice of the weights Wi is necessary From general consid¬ 
erations 0 only two independent boundary 0(a) counterterms involving the gauge helds 
/ d^xFofc(x)Fofc(x) and / d^xFki{x)Fki{x) at Xq = 0 and Xq = T are expected. Since we 

^For improvement of some physical quantities arising from local Ward identities a particular choice of 
Wi may not be necessary. 
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now have many types of loops at onr disposal there are various ways to chose the weights 
to achieve our goal. A rather natural way, which we call Choice A, is to adjust the 0(a) 
counterterms just through the plaquettes: 

( Csigi) if ii®® completely on one of the boundaries, 

Ct{gQ) if C just touches one of the boundaries, (2-13) 

co(fi'o) otherwise. 


and for the other classes 5*, z = 1, 2, 3 we simply set: 


Wt{C,gl) = { 

0 if C lies completely on one of the boundaries, 

Ci(5'o) otherwise. 

(2.14) 

With this particular choice for complete 0(a) improvement we hnd to tree level 



Cs(0) = 

(2.15) 


Ct(0) = Co(0) -|- 2ci(0), 

(2.16) 


independently of the (smooth) boundary conditions. The necessity for the extra term in 
Ct(0) involving the coefficient of the rectangles is explained in Appendix A. 

From the discussion in Appendix A it will also be clear that another admissible special 
way to set things up, which we call Choice B, is to keep eg. ( |2.13|) but change eg. (|2.14|) to 

{ 0 if C lies completely on one of the boundaries, 

|cj( 5 fg) if z = 1 and C has exactly 2 links on a boundary, (2.17) 
Ci{go) otherwise. 


Then eg.(|2.15D remains valid but instead of eg.(|2.1ti|) we must set 


Ct(0) = Co(0). 


(2.18) 


In other words in this case one of the 0(a) boundary counterterms at tree level is adjusted 
by the coefficient of the rectangular loops of space-length 2 lattice units. The reason for 
mentioning this choice will become manifest below. 

We must still specify the boundary conditions for the gauge held. As in the determi¬ 
nation of the running coupling 0 we chose 

Uk{x)\xo=o = , Uk{x)\^g=T = , (2.19) 


7 







where C, C are diagonal traceless antihermitian matrices. As explained in ref. the pres¬ 
ence of non-trivial boundary conditions induces a constant chromo-electric background 
held. In the continuum limit the background held is given (up to gauge transformations) 
by 


Aq{x)= 0, AA.(a;) = 6 (a;o), fc = 1,2,3, 


h{xo) = h {T - xo)C XqC 


On the lattice one will have a solution to the held equations of the form 


( 2 . 20 ) 

( 2 . 21 ) 


Uo{x) = 1, Uk{x) = 'O(xo), k = 1,2,3, 


( 2 . 22 ) 


but for actions having ci 7 ^ 0 and generic choices of (admissible) boundary weights it is 
not possible to solve the lattice equations of motion analytically and one must be content 
with numerical solutions. Although this is of no relevance for numerical simulations, for 
perturbative computations it is more practical to have analytic solutions at tree level. 
This is where the special Choice B comes in to play since it is easy to check that in this 
case 

l/(a:o) = exp6(a:o), (2.23) 

with b as in eq. (|2.21|) , is an exact solution (in the limit (70 = 0). On the other hand for 
Choice A, although the above form of C is a solution to the equations obtained by the 
variation of links well inside the cylinder, it fails to be a solution for equations obtained 
by the variation of spatial links at time 1 (or at time T — 1). Thus for non-vanishing 
Abelian background held it is clearly advantageous to use Choice B for the analytical 
computations. 

The Feynman rules are now easily derived by expanding around the background held 
solution 

Uo{x) = exp[goqo{x)], Uk{x) = V{xo) exp[goqk{x)]. (2.24) 


For the validity of the perturbative expansion it is necessary to know whether the back¬ 
ground held is the unique absolute minimum of the action (up to gauge transformations). 
For the case of the Wilson plaquette action this was proven [Q provided the diagonal 
elements of C lie in a certain so called fundamental domain. The extension of the proof 
of this theorem to the general class of actions considered here has to our knowledge not 
been achieved. Klassen [|^ has however investigated the question numerically and found 
no evidence to the contrary for a large class of actions. In fact the only cases where 
Klassen found that the global minimum was not given by eqs.(|2.23|) and (p.21|) were those 









with coefficients not in the positivity domain P alluded to above. We hope therefore that 
perturbation theory remains trustworthy. As a consistency check in our computations we 
monitored the eigenvalues of the quadratic gauge operator (the inverse free propagator) 
and checked (for various covariant gauges) that it did not develop a negative eigenvalue. 
This held even for the action RG3 which is not manifestly positive. 

Concerning gauge hxing we used the covariant gauge hxing specihed in ref. 0 


S^ = Xo{d*q,d*q) 


(2.25) 


with R 

( for 0 < Xo < T, 

{d*q)af){x) = I L -3 y)apda0 for Xo = 0 , (2.26) 

I 0 for xq = T. 

The inverse gauge propagator for the case of zero background held is given in Appendix 
B. Note that the simple choice of gauge hxing given above has the disadvantage that there 
is for generic q no choice of the gauge parameter where the propagator becomes diagonal 
with respect to the Lorentz indices. It may be possible to introduce other choices of gauge 
hxing where this is the case but we did not investigate this possibility further. In all cases 
we obtained the free propagator numerically although analytical expressions are probably 
available for the case of vanishing background held. The expressions for the vertices are 
rather lengthy and we do not give them here, they can be obtained from the authors upon 
request. 


2.3 The fermion action 

As mentioned above for fermion action we consider only the on-shell 0(a) improved Wilson 
action which is obtained by adding the Sheikholeslami-Wohlert interaction term (with 
coefficient Csw) H]- That is in eq. o Sw is the standard Wilson fermion action with 
r = 1, 

S'w[C,^,^] = '^'4}{x){D + mo)'ilj{x) + S'wb, (2.27) 


where mo is the bare quark mass and 

o=5i:n4v;+v„)-v;v„}, 


(2.28) 


^The reader is requested to consult ref. ^ for the specific notations. A more symmetric way to fix 
the gauge is discussed in the paper by Narayanan and Wolff M. 
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is the lattice Wilson Dirac operator. V* are the forward and backward lattice covari¬ 
ant derivatives 


= X^U^{x)'4){x + fi)- ipix), 

Vl'4){x)= '4){x) - X^U^ix - (i)^{x - fi). 
The additional U{1) constant held 

1 if /i = 0 , 




3*®/^ if /i > 0 , 


(2.29) 

(2.30) 

(2.31) 


corresponds to a free phase in the spatial boundary conditions of the fermions. Finally 
S'wb in ( p.27| ) is an additional boundary term []^ which is needed to remove 0(a) lattice 
artifacts in the SF framework: 


^wb = {ct-l)T,x{^xo,ii’{x) 'ipix)-Uo{x-6) ^P+'ipix-O) 
+(5j;o,T-i^(a;) i^{x) - Uo{x)P-ip{x - 1 - 0) 


(2.32) 


where P± = 1(1 ± 70 ). 

5sw is the Sheikholeslami-Wohlert [H term: 


SswP.ip.ip] = i XI 'ip{x)a^^P^^{x)il){x), 


(2.33) 




where the held is given explicitly by 

= |{[ U^{x)U,{x + fi)Uf.{x + OyU^ixV 

+ U^{x)Uf,{x + 9 - - iiyUf,{x - p,) 

+ U^{x — ppUy{x — p — 9pU^{x — 9 — p)Uy{x — 9) 

+ f/^(a; - 9pUf,{x - 9)U^{x + p - 9)Uf,{xp] 

[/x^z/]}. 


(2.34) 


The boundary conditions for the fermion helds have hrst been specihed in ref. [|I 8 


: m 


all our computations we have set the quark held on the boundaries to zero. 
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A general introduction to on-shell 0(a) improvement of lattice QCD with Wilson 

In particular there it has been discussed how the 


quarks can be found in refs. [21, 22 


renormalization procedure can be carried out in a way consistent with 0(a) improvement. 
In the following we adopt the conventions and notations as introduced in this reference. 


2.4 Composite operators of dimension 3 

On-shell 0(a) improvement requires also the introduction of 0(a) counterterms for the 
composite helds. Here we restrict attention to a few gauge invariant composite operators 
which are bilinear in the quark helds. Assuming Nf > 2 mass degenerate quark havors 
we consider the local isovector helds 


Kix) = ^(a;)7^ir“V’(a:), 

(2.35) 

= ^(a^)7/.75|r“^/>(a;), 

(2.36) 

P“(x) = 

(2.37) 


(2.38) 

Using these dehnitions and following ref. [Q the corresponding improved vector, axial 
vector and tensor currents may be parametrized as follows. 


(2.39) 

(Ai)l= /1J + ca(9o)AC“, 

(2.40) 

+ cr(gl)a{a^v‘ - 

(2.41) 

where d denotes the symmetric lattice derivative d^f{x) = ^[f{x + afi) — f{x — afi)]. 

Throughout the paper we will use a mass-independent renormalization scheme which 
is compatible with 0(a) improvement. In such a scheme the renormalized coupling and 
quark mass are given by 0 

9l= 9oZg{glafi), 

(2.42) 

ruR = mqZ^{g^,an), 

(2.43) 


® The Dirac gamma matrix conventions are as in Appendix A of ref. ||^ . 
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where /i is the renormalization scale. The parameters go and rriq are dehned as follows 
9o= + hi9i)a'mq}, (2.44) 

fhq = mq[l + b^{g^)amq], mq = mo - me, (2-45) 


where rric is the critical bare quark mass. The critical hopping parameter is thus given 
by i^c = 1/(8 + 2amc) . 

With these dehnitions the renormalized improved bilinear operators take the form 
(X = V,A,P,T) 

Ar = Zx{gl, a/i) [1 + bx{gl)amq\Xi. (2.46) 


For the Wilson plaquette action a non-perturbative determination in the quenched 
approximation of Z^, by can be found in |^; of Ca in [Q; and of Cy in . 


2.5 Computation of the improvement coefficients 

One hrst considers the massless theory. In this case on-shell correlation functions involving 
only the composite helds are 0(a) improved provided the improvement coefficients Csw, 
Ca, Cy and ct are properly chosen as functions of the bare coupling go. In the Schrbdinger 
functional framework it is, however, usually easier to perturbatively compute correlation 
functions also involving the near-boundary quark helds ( and 0, 0._O(a) 

improved correlation functions are then obtained with the renormalized helds Cr, Cr cind 
Cr5 Cri which are all related to the bare helds by the same renormalization factor, e.g. 


Cr = Z<^{gl,ag.)[l + 6^(^o)amq]C. 


(2.47) 


The particular set of correlation functions which we have chosen to compute was 

Our 


precisely as for the case of the Wilson action described in detail in refs. |]12|, ^ 


calculation determines of course only the 1-loop correction to the tree level values of 
the improvement coefficients, i.e. the quantities in the perturbative expansion for a 
general improvement coefficient cx'- 


( 0 ) , ( 1 ) 2 , ( 2 ) 4 , 

Cx — Cx + Cx (70 + '^x 5'o + • • • 


(2.48) 


Only for the computation of is it necessary to have a non-vanishing background 
held. Once one knows this coefficient the others can be determined by computing corre¬ 
lation functions in the absence of a background held, for which the programs are consid¬ 
erably less CPU time consuming. In the case of vanishing quark mass (mRL=0) we have 
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collected data at T = 2L, for 6 = 0 and 6 = 1, with vanishing backgronnd held and at 
T = L, only for 6 = 0, with non-vanishing backgronnd held. 

After calcnlating the improvement coefficients which are necessary in the massless 
theory one can proceed to the determination of the 1-loop contribntion to the various 
6 -coefficients that were introduced above and are required for the theory with massive 
quarks. These coefficients have an expansion in powers of Qq of the same form as the 


expansion (|2.48|) for the c-coefficients. In the massive case one has to study quark correla¬ 
tion functions at hxed values of m-^L. The set of correlation functions that we considered 
is exactly the same as described in detail in ref. |@. Also in this case the computation 
can be done at vanishing gauge background held. Data were collected for two values of 
turL, namely 0.1 and 0.5 and 6* = 0, by taking T = 2L. 

We remark that we did not directly determine the coefficient for the isovector 
scalar operator S°‘ since it can be determined by the relation 


k' = - 26 «. 


(2.49) 


In the quenched approximation one can prove the non-perturbative relation 6s = ~26in 
^ from which the 1-loop relation above follows for the fully interacting theory. We also 
observe that the improvement coefficient 6g (see Eq.(2.44)) vanishes at tree level and at 
one loop order it is independent of the pure gauge action: 6^^^ = 0.012000(2)W p7| . 

It should also be noticed that the values of and 6^^^ actually needed for the 0(a) 
improvement of correlation functions involving boundary quark helds depends in general 
on the precise form of the pure gauge action. This involves not only the diherent choices 
of the Wilson loop coefficients c, {i = 0,1,2,3), but also the particular choice made for 
the SF setup. Our results for and 6^^^ refer specihcally to the above mentioned Choice 
A for all the considered pure gauge actions. 


3 Analysis Method 

We describe in this section the method that we adopted to extrapolate the asymptotic val¬ 
ues of the various improvement coefficients from the 1-loop results obtained on a sequence 
of hnite lattices. 

The 1-loop results were produced by Fortran codes in double precision. We have 
checked the gauge invariance of the results produced with our codes by adopting different 

® A brief outline of the proof, which is due to M. Liischer, can be found in the last section of ref, . 
Please note that in this paper the factor on the rhs of eq.(113) should be (1 -I- 2bmamq)~^. 
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covariant gauge fixings within the family of eq.( p.25|) . Checks were also done against pre¬ 
existing safe codes in some limiting cases, like the case of the standard Wilson plaquette 
action or the case of vanishing background field. Further checks about the correctness of 
the codes and of the whole theoretical setup were obtained a posteriori by the observed 
consistency of independent determinations of the same improvement coefficients. 

The effects of rounding errors and their dependence on the lattice size have been 
investigated by comparing, for a limited set of cases, results obtained by running the 
same code in single, double and extended precision. From this investigation we have 
derived a formula which gives a rough estimate of the rounding-error effects as a function 
of the lattice size. The estimated rounding errors were conservatively propagated through 
the whole analysis procedure described below. In all the cases, when the analysis output 
was clearly affected by rounding errors, our estimate was found to be realistic or even 
conservative. 

The starting point of the analysis procedure are the (gauge invariant) correlation 
functions /x(T; T/L; 6 *; turL), evaluated on a sequence of lattices of volume x T at a 
given fixed value of ttirL and with the spatial boundary conditions for the fermion held 
specihed by 6, as explained in section 2. Here X labels the different quark correlation 
functions within the set considered in refs. HT^ ^6 


Taking suitable linear combinations of these correlation functions, possibly involving 
different values of 6 and m-^L, enables us to construct lattice approximants of the improve¬ 
ment coefficients or some simple linear combinations of them, which are the quantities 
that we will consider from now on. The precise form of these quantities is not given here, 
because, as already mentioned, they are exactly the same as specihed in |]^, for the 
determination of K^\ and and in for all the remaining improvement 

coefficients introduced above. 

Each sequence {Q} of lattice approximants has, by construction, a well dehnite limit 
as L ^ oo and a large L expansion of the form: 


Q(L; TjL) = Qo + ElCi + 


(3.1) 


i=l 


where the lattice artifacts coefficients Qi and g* depend on the value of T/L. We are 
of course interested in the asymptotic value (Qo) of the sequence, which is obtained by 
suitable extrapolation to the limit L —> oo, performed at a hxed the value of T/L. 

In order to get the asymptotic value of the various sequences of interest we have used 
the method of |2^, suitably adapted to the case of functions of L of the form (|3.1|) . The 
main point of the method is to perform some consecutive blocking steps on the initial data. 
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represented by the sequence {Q}. In the hrst step, for each value of L, one constructs the 
linear combination of the data at L and the two contiguous values of the lattice size {L — 1 
and L + 1, in the case T = 2L) that does not contain terms of order L~^ and log(L)L“^. 
In the second step, the terms of order L~^ and log(L)L“^ are analogously removed and so 
on. At each step the number of data in the resulting sequence is of course decreased by 
two units and the rounding-error effects are enhanced, which limits the optimal number 
of steps to be performed in practice. This method allows for a safe, although often 
conservative, estimate of the systematic error on the asymptotic value of the sequence 
under study. To this aim it is important to properly propagate the estimated rounding 
errors on /x(T; T/L; 0; mRL) when constructing all the needed linear combinations. Our 
practical use of the method closely followed the original proposal and detailed description 
in [^1 


apart from some details which are discussed below. 


The data for the correlation functions /x(T; T/L; 0; mj^L) were collected for values of 
L ranging from 4 to L^ax, as indicated in table |^. Because our data have been obtained 
using only double precision and with a value of L^ax = 32, the uncertainties on the 
values of improvement coefficients that we quote for LW, RGl, RG2 and RG3 actions are 
systematically larger than the uncertainties obtained in |^, for the standard plaquette 


action. In addition, we had also to adopt a slightly modihed blocking-£t procedure, in 
some cases where a straightforward application of the above mentioned blocking procedure 
turned out to lead to a large (order 10%) estimated value of the relative error on the 
asymptotic quantities we are interested in. 

We realized, indeed, that when the lattice artifacts of a given order j in the asymptotic 


expansion (Rl|) , i.e. 

[Qj + Qj log(T)]L"-’ 

are numerically very small (in comparison with higher order lattice artifacts) for all the 
considered values of L, removing them through a suitable linear combination strongly 
enhances rounding-error effects while hardly changing the estimate of the asymptotic 
value {Qo in (|3.1| )). In such a case we found it convenient for the purpose of improving 
the precision on Qo, to resort to the following hybrid blocking-fit procedure: 


perform standard blocking steps for terms of order from 0 to j — 1; 


• using suitable linear combinations of the reduced data sequence at step j — 1, de¬ 
termine separately the values of Qj and qj, say Qj and qj, with some (in practice 
rather large) relative uncertainty 6Qj and 6qj; 

• subtract from the reduced data sequence at step j — 1 the quantity [Qj+qj \og{L)]L~^ 
and perform standard blocking steps for terms of order j + 1, j + 2, ..., stopping as 
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usual when the number of data in the sequence gets too small or the rounding-error 
effects too large 


• repeat the procedure of the latter item by replacing Qj and qj with Qj ± 6Qj and 
Qj ± 6qj, respectively, taking mean and extreme values in a way that is equivalent 
to the error estimate adopted for standard blocking steps. In this way one gets a 
conservative estimate of the uncertainty induced by SQj and Sqj. 


This blocking-£t procedure is equivalent to the simple blocking procedure of in the 
ideal case of inhnite precision data, but turns out to be more convenient when rounding- 
error effects are not negligible and lattice artifacts of a given order j are very small. 

We remark that requiring 0(a) improvement of all the correlation functions introduced 


in ref.s |^, ^ provides a number of relations between the various improvement coefficients 
larger than the number of the improvement coefficients themselves, allowing then for 
several consistency checks. Further consistency checks were done by comparing the values 
of the b-coefficients at = 0.1 and turL = 0.5. 

After checking the consistency among the asymptotic values obtained from different 
lattice approximants of the same improvement coefficients, we have sometimes made use 
of the fact that a particular linear combination of these approximants showed small lattice 
artifacts, leading to a more precise determination of the asymptotic value than the single 
approximants themselves. 


4 Results and Conclusion 


In table ^ we give our hnal result obtained by the method in the previous section. This 
table contains a synthesis of all our 1-loop results, including a comparison with the case 
of the standard plaquette gauge action (label “Plaq.”)|^, 

We hrst notice that tadpole contributions to the denoted cl^ in the table, give 
about 90% of the complete 1-loop contributions for all actions considered here. Therefore 
the value of Cgw taken by the CP-PACS collaboration for their full QCD simulation with 
RGl gauge action [Q is very close to the full one-loop value up to order gQ, and the order 
(7o difference can be calculated from the table: + O.OOS^fQ -|- 0 (( 7 g) . 

Although only 3 choices of RG improved gauge actions are considered, it seems that 
RG improved gauge actions generally give a coefficient which is a factor 2 to 2.5 smaller 
than that of the plaquette action, while for the perturbative improved action (LW) the 
reduction factor is only about 1.35. This tendency has already been found in the finite 
part of 1-loop renormalization factors for various quantities [^. Other improvement 
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coefficients, the Cx^’s and &x^’s, as well as and c[^'^ for the setup denoted as Choice 
A, are also given in table |I[ The same tendency for the size of 1-loop coefficients is also 
observed. 

It is clear, however, that the smallness of 1-loop coefficient for the RG improved action 
does not imply the smallness of the lattice artifacts for the action. To get some impression 
of these we have to evaluate lattice artifacts directly from physical observables for the 
corresponding action. For that purpose, following ref. |]^|, we consider the unrenormalized 
current quark mass m, dehned through the PCAC relation (see eq.(6.13) of ref. [0 ), which 
is expected to vanish up to terms of order a? ai K = K^. If we expand m such that 


am = ro + viC^gl + 0 ( 5 ( 3 ) 


(4.1) 


at K = Kc and Xq = T/2, the value of ri represents the magnitude of the remaining cutoff 
effect for the particular gauge action at 1-loop order, since the tree-level contribution tq 
is independent of the pure gauge action. In £g. ^ we have plotted ri for various actions 
as a function of a/L for T = 2L and 6 = 0 and vanishing boundary gauge field (in which 
case the tree-level contribution actually vanishes). At large values of a/L it is observed 
that this 1-loop lattice artifact is indeed smaller for the RG improved gauge actions than 
for the LW action, for which it is still smaller than for the plaquette action. We remark 
however that one cannot draw general conclusions from this simple exercise; the pattern 
may be different for other correlation functions and, needless to say, at a non-perturbative 
level we have no statement to make at all. 

In the last box of the table the 1-loop value of the critical hopping parameter 
dehned by + ■■■■, is also given and compared with the results from 

perturbation theory on an inhnite lattice . They agree with each other within 5 digits, 
and this fact also supports the correctness of our calculations. 

Recently the coefficients cx and 6x have been calculated by Taniguchi and Ukawa 
im using a completely different method. The two sets of results agree well, although 
our quoted errors are generally larger than theirs. The good agreement between the 
central values of the two methods indicates, however, that our error estimates are rather 
conservative. 
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Appendix A 


In this appendix we explain the allowable weights of the contributions to the gauge action 
near the boundary in our setup to avoid order a effects. We remind the reader that in 
this paper we are not demanding the absence of O(a^) effects. As explained in ref. 0 
there are only two possible boundary counterterms at order a in the pure gauge theory, 
which are obtained by summing any local lattice expressions for the helds 


tr 


l^FokFok'^ 


and a^tr 


{^FijFij 


(A.2) 


over the xq = 0 and xq = T hyperplanes. 

The hrst step is to ensure that the approach to the classical continuum limit is O(a^). 
Following the discussion of ref. [^| we recall that any given smooth continuum gauge held 
can be arbitrarily well approximated by lattice gauge helds U^{x) 


U^{x) = T expja f dtA^{x + [1 - t]a/i)} • 


(A.3) 


One can then derive asymptotic expansions for the lattice operators. For example for 
C{Ci) = ReTr[/ — U{Ci)] with Ci G Si one easily hnds for the hat loops in the /i — z/ plane 


C{Ci) = -kitt^ ti F^^{xc}F^^{xc) + 0(a®), i = 0,1 


(A.4) 


where the constant fcj depends on the class 

ko = ,ki = 2. (A.5) 

and Xc is the “center of the loop” (otherwise one obtains in general O(a^) appearing on 
the rhs of Eg. ( A-41) above). Obviously there are similar expressions for the non-planar 
loops. 

The approximation of the integrals by sums in spatial directions can incur only O(a^) 
errors because of the periodic boundary conditions. However approximations of the in¬ 
tegral by sums in the time direction may result in additional ehects of order a through 
boundary terms. For example for a smooth function f{z) dehned on the region [0, T] and 
dividing the region in M = T/a segments, we have 

M rp 

a'^f{ma—\a) = / dzf{z) + 0{a^), (A.6) 

m=l F 


M-l 

m=l 


£dzf(z)-^{m + f{T))+0{a^). 


(A.7) 
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For planar loops extending in the time direction for which the time coordinate of the 
center Xc takes half integer values we can invoke Eq.( A. 6 ) to see that their contribution 
to the action will be as in the continuum up to terms of order a^. However for those loops 
where Xc takes integer values, which is the case only for the rectangles extending over two 
time units, we see from (|A.7|) that the contribution to the action from these rectangles 
will have an 0(a) dehcit with respect to the classical continuum expression. For instance 
at the boundary T = 0 we are missing a contribution —( 2 ci/( 7 q) / dxtr Foa;( 0, x)Fofc(0, x) 
which can be compensated by adding extra terms involving plaquettes and/or rectangles 
extending only one unit in the time direction. Two admissible choices are the special 
setups A, B specihed in subsection 2.2. 

Eq. (|A.71) also explains the fact that Cs(0) = but this is not relevant for our particular 
boundary conditions for which Fjk — 0 . 


Appendix B 


In this appendix we give the form of the free inverse gauge propagator for non-zero ci and 
C 3 in the case of vanishing background held. Introducing the Fourier transformed helds 
through 


Qoix) 

= L ^'^F^^qo{p,Xo), 

p 

(B. 8 ) 

Qkix) 


(B.9) 


p 


the quadratic part of the gauge action including the gauge-hxing term takes the form 

T-l 

>^£uge = Z! H t^Q,i{-P,Xo)q^{p,yo)K^^{p;xo,yo). (B.IO) 

P xo,yo=0 


Explicitly 


A'oo(p; Xo, yo) = Sxo,yo [(co + 6ci 8C3)p^ - Cip^ -h C3(p^ - (p^)^) 

+ •^o{2 — <^a:o,o(l ~ ^p,o) ~ ^xo,r-l} 

+ (<^3;o-l,yo F 5x0+1, I/O )(cip2 - Ao) 

T cidxQ^yQ(^xo,T—i T 5xo,o)(p P 

Fko{p] xq, yo) = ipk ( 5 x 0 ,j/o ~ 5xo-i,j/o){co + 5ci -f 8 C 3 — Aq — ciPk^ + C3{pk^ 


(B.ll) 

P^)} 
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Kokiv'^xo^yo) 

Kki{v\XQ,yo) 


+ Ci((5j;q + 1 ,J/q 5x0-2,J/o) 

-|- Cl ( 5 j;q ^xo — l,yo^xo,l )(1 Pk 5 fec) 

= —Kko{p;yo,xo), 


~ ^XQ,yQ^kl 
2 


(B.12) 

(B.13) 


|^(co + 8 ci + 4 c3)p + 2 co + lOci + I 6 C 3 

ci(2pfc^ + + pVfc^) + C3(2pfc^ + P^ - (p^)^ + pVfc^) 

Uo - Co - 8Ci - 6C3 + Ci(pfc^ + pi^) + C3(p^ - pk^ - pi" 

2 , „ 2 oa2^ 


^XQ^y^PuPl [ 

(5x0+1,j/o 5" 5x0— 1 , 2/0 ) |5a;/{co + 4ci + 8 C 3 — ciPk^ + C2,{pk^ — 2p^)} 


+ 

+ 


CzPkPl 


(5x0+2,2/0 4" ^xo-2,yo)SklCl 


Cl<5kl^xo,yo(^xo,T—l 4“ 5 ^ 0 , 1 ) (1 Pk ^fec), 


(B.14) 


where 


Pk = 2sinpfc/2, 


P 



A:=l 


P 


EPk^ 


k=l 


(B,15) 


For the boundary coefficient d^c appearing above, we have d^c = 0 for Choice A and 
dbc = 1 for Choice B. 
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vs. a/L 


0=0 



a/L 


Figure 2: Remaining cutoff effect at 1-loop in the PCAC mass at xq = T/2 and K = Kc 
for various actions, for the case T = 2L, 6 = 0 and vanishing background held. 
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Table 1: Improvement coefficients at 1-loop for various gluon actions. The results for the 


standard plaquette action are taken from refs. |]T2 
with T = 2L except d 


^^ ^_[] . All results are obtained from data 

which is calculated from T = L data. Tadpole contributions for 
are also listed. The critical hopping parameter (above) is also compared with 
the result (below) obtained from the usual perturbation theory in inhnite volume [p0|. 


gauge action 

Plaq. 

LW 

RGl 

RG2 

RG3 

Cl 

0.0 

-1/12 

-0.331 

-0.27 

-0.252 

C3 

0.0 

0.0 

0.0 

-0.04 

-0.17 



Lftiax 

32 

24 

30 

24 

24 

Ai) 

'^SW 

0.267(1) 

0.196(6) 

0.113(3) 

0.119(5) 

0.109(5) 

^sw 

0.25 

0.183 

0.105 

0.110 

0.096 

Lynax 

48 

32 

32 

32 

32 




xCf 



Ji) 

Ca 

-0.005680(2) 

-0.004525(25) 

-0.002846(11) 

-0.003017(12) 

-0.002805(20) 

Ji) 

Cy 

-0.01225(1) 

-0.0103(3) 

-0.00730(20) 

-0.00757(26) 

-0.00709(20) 


0.00896(1) 

0.00743(7) 

0.00505(10) 

0.00526(15) 

0.00496(12) 

!>'£’ 

-0.07217(2) 

-0.0576(11) 

-0.0382(8) 

-0.0395(15) 

-0.0353(12) 

uW 

0.11414(4) 

0.0881(13) 

0.0550(4) 

0.0572(6) 

0.0510(5) 

Oy 

0.11492(4) 

0.0884(26) 

0.0551(19) 

0.0574(19) 

0.0510(21) 

Op 

0.11484(2) 

0.0889(14) 

0.0558(9) 

0.0584(10) 

0.0528(8) 

l(1) 

0.14434(4) 

0.1152(22) 

0.0764(16) 

0.0790(30) 

0.0706(24) 


0.10434(4) 

0.0790(25) 

0.0477(12) 

0.0502(19) 

0.0444(15) 


-0.06738(4) 

-0.0505(23) 

-0.0243(11) 

-0.0275(18) 

-0.0252(11) 

cl 

-0.01346(1) 

-0.0122(4) 

-0.00661(21) 

-0.00612(23) 

-0.00186(20) 

Lynax 

32 

32 

32 

32 

32 




xCpX 10-3 



Ki}) 

6.329891(3) 

4.716253(1) 

2.760894(1) 

2.911363(2) 

2.593490(1) 

KW 

6.3300(1) 

4.7163(1) 

2.7609(1) 

2.9114(1) 

2.5935(1) 
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